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1. Motivation and Objective 

True dissipation occurs mainly at the highest wavenumbers where the eddy sizes 
are comparatively small. These high wave numbers receive their energy through 
the spectral cascade of energy starting with the largest eddies spilling energy m o 
the smaller eddies, passing through each wavenumber until it is dissipated at the 
microscopic scale. However, a small percentage of the energy does not spi con in 
uously through the cascade but is instantly passed to the higher wave numbers . 
Consequently, the smallest eddies receive a certain amount of energy almost imme- 
diately. As the spectral energy cascade continues, the highest wave number needs 
a certain time to receive all the energy which has been transferred from the largest 
eddies. As such, there is a time delay, of the order r, between the generation of 
energy by the largest eddies and the eventual dissipation of this energy. 

For equilibrium turbulence at high Reynolds numbers, there is a wide range where 
energy is neither produced by the large eddies nor dissipated by viscosity, but is 
conserved and passed from wavenumber to higher wavenumber. The rate at which 
energy cascades from one wavenumber to another is proportional to the energy 
contained within that wave number. This rate is constant and has been used in 
the past as a dissipation rate of turbulent kinetic energy. However, this is true 
only in steady, equilibrium turbulence. Most dissipation models contend that t e 
production of dissipation is proportional the production of energy and that the 
destruction of dissipation is proportional to the destruction of energy. In essence, 
these models state that the change in the dissipation rate is proportional to the 
change in the kinetic energy. This assumption is obviously incorrect for the case 
where there is no production of turbulent energy, yet energy continues to cascade 
from large to small eddies. If the time lag between the onset on the energy cascade 
to the destruction of energy at the microscale can be modeled, then there will be a 
better representation of the dissipation process. Development of an energy cascade 
time scale equation will be discussed in Section 2.2. 


2. Work Accomplished 
2.1 Mean Flow Equations 

For incompressible flow, the equations for continuity of mass and axial momentum 
are written as aTT 
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where, uiuj is the turbulent Reynolds stress tensor. Using the' eddy viscosity 
concept 1 2 , the Reynolds stress may be related to the mean strain rate and a turbulent 
eddy viscosity, 

-m 

The turbulent viscosity given in the above equation can be interpreted as a measure 
of the turbulent kinetic energy existing in the flow times a local length scale, 
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v t = ck^l 


where c is an arbitrary constant. The definition used for this length scale is the pri- 
mary discriminating factor between turbulence models and determines the number 
of equations which need to be solved. The length scale can be written in terms of 
the turbulent kinetic energy, fc, and its dissipation rate,c, 
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Now, the momentum equation can be written as 



where, 
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and c p is a constant. The k and e are determined by solving transport equations 
for the turbulent kinetic energy and the dissipation. These equations and the mod- 
ification to the “standard” dissipation equation will be discussed in Section 2.2. 


2.2 Turbulence Equations 

Taking the Navier-Stokes equations and multiplying by w, = Ui + u, and then 
taking the time average yields an equation that contains both the energy equation 
for the mean velocity, Ui, and a mean energy equation for the fluctuating component 
of velocity, u,-. Eliminating the energy equation for the mean velocity yields the 
turbulent energy budget in the absence of a pressure gradient field 3 , 
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The last term in equation (3) is the rate at which viscous forces perform deformation 
work to the fluctuating strain rate, defined as 
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This term is the true viscous dissipation rate which drains kinetic energy from the 
system, e = This dissipation rate is important to the overall turbulence 

structure. The next to the last term in equation (3) serves to transport kinetic 
energy between the mean flow and the turbulence. It is generally referred to as the 
kinetic energy production term. Writing the mean strain, Sij , as 
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the production term for incompressible flow is defined as 
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The last set of terms to be modeled represents the energy redistribution by viscous 
forces and are contained within the parenthesis on the right hand side of equation 
(3). From the definition of the fluctuating strain rate, 
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and using a mean gradient hypothesis, the triple correlation can be rewritten as 
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where a k is a constant of the order one. Combining all these modeling assumptions 
yields the high Reynolds number form of the turbulent kinetic equation, 
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of the following form has been widely used for modeling of the 



The production term, P t has been assumed to be proportional to the turbulent 
kinetic energy production term, equation (4). This approximation places a di- 
rect correlation between the production of turbulent energy and the production of 
dissipation. Although this may be true at equilibrium, where by definition the dis- 
sipation equals the turbulent energy production, there are some deficiencies in this 
statement for the more general case. According to Lumley 1 , the modeled dissipation 
rate equation should be written as, 
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The production term, P e , now represents the rate at which energy arrives at the dis- 
sipative wavenumbers through the energy cascade starting with low wavenumbers. 
In equation (5), 1/S represents the characteristic time for energy to be transported 
to the dissipative wave numbers. 

In the following equation, the term on the left hand side and the first term on the 
right represent the transport of the inverse time scale by advection and diffusion, 
respectively. The remaining terms are the production from the mean straining 
forces within the large scale mean flow structures and the dissipative losses at the 
microscale, 
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The production and loss terms have been scaled by a characteristic time, T, which 
in this study has been taken to be a weighted function of the time scale for the large 
eddies, r, 
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Since most energy enters the cascade at scales comparable to the integral scale, l, 
the characteristic time scale of the large eddies dictates the rate of energy transfer 
and is defined as r = £ for nearly isotropic turbulence. Using the approximation 

for the integral length scale as, l = the time scale becomes 


r 



which is valid even in nonequilibrium flows. Letting Cb be a function of mean and 
turbulent quantities, the inverse time scale equation can be written as 


DS 

Dt 





where , 

Cb — i + ^y/SjjSjj 

Determination of the modeling constants is discussed in the next section. 

2.3 Model Constants 

In this paper, only high Reynolds number turbulence has been considered. The 
model constants are examined for the grid turbulence case where the mean velocity 
gradients and strain rates vanish. For this case, turbulent intensity decays with the 
following proportionality, 
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so that the dissipation rate must decay as 
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Here, n is the decay rate and is typically in the range 1.1 <n< 1.3. Substituting 
these relationships into the inverse time scale equation yields 
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After integration, this equation is 
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where a is a constant of integration. For there to be an inverse scaling relationship 
between S and t, the exponent in the above equation must reduce to -1 for the gnd 
turbulence case. Since the coefficient C B reduces to 1.0 in grid turbulence, then 

cb — n. 

In this study, n has been set to 1.1 and held constant 4 , which sets the value of c B . 
Now using the fact that, initially, 

So ‘° = in 

the coefficients c A and c D are related as follows 
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The coefficient has been evaluated for equilibrium turbulence where, where 
P k = e. Rewriting equation (2) as 
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and using equations (1) and (4), this coefficients is 
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After the above analysis, the only adjustable parameters are and the a coeffi- 
cients, which should be of the order one. The complete set of optimized coefficients 

is given in Table 1. 
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Table 1. k-e-S Model Coefficients 
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2.4 Free Shear Flows 

The three turbulent transport equations are solved implicitly using the finite 
difference technique of Spalding 5 for parabolic equation systems. An initial profile 
is imposed upon the first marching plane for the turbulent kinetic energy and the 
dissipation, and S is initially set to S, = In the mixing layer and jet cases, 
the initial plane was equally divided into a uniform velocity field and still air. For 
the case of the plane wake, this initial profile came from a flat plate calculation 
computed with the same solution technique. Starting with the initial plane, the flow 
field is developed by marching the transport equations downstream and applying 
the following boundary conditions, 

Planar Jet 

Centerline (y = 0) f* = 0 > = 0 

~ _ S o liC&fcCfy = 0 + ^((*'+* , t)$£) 

O J/=0 ~ c h kU v=0 +eAx 

Far Field (y = y m ax) A: = 0,e = 0,<S = 0 
Round Jet 

Centerline (r = 0) §£ = 0 , = 0 

^*•=0 — c t kU r , o+eAx 

Far Field (r = r max ) fc = 0,e = 0,S = 0 
Planar Mixing Layer 

Centerline ( y = y m in) m ' " T = 0-015, — = 3.5 x 10“ 6 

w y — ymxn y — ymin 

C — 1 

°y=ymm — 

Far Field (y = y max ) k = 0,e = 0,S = 0 
Planar Wake 

Centerline (y = 0)fj=0,|j=0 

SoliCi,kU v= o + -&; (( l/ + ,/ «)^) 

«>y=0 — C k Jfct/ v = 0 +<Al 

Far Field (y — y max ) fc = 0,e = 0,5 = 0 

The solutions have been checked after several hundred steps to insure that the 
profiles have become self-similar. 

2.4.1 Planar Jet 

The three equation k-e- S model performed very well for planar jet flow when 
compared to the experimental data of Gutman and Wygnanski 6 , Bradbury 7 and 
Heskestad 8 . Mean velocities profiles given in Figure 1 are well predicted by both the 
k-e- S model and the standard k-e model, which is reassuring since the coefficients 
for both of these models were optimized for the planar jet case. The spreading rates 
predicted by both models are within the experimental range and are given in Table 
2 . 

There is a slight discrepancy between the k-e- S model and the standard k-e 
model for centerline kinetic energy as can be seen in Figure 2. However, both 
predictions are well within the range of experimental measurements, which show 
a considerable spread near the centerline of the jet. Both predictions are closely 
aligned with the data of Bradbury 7 . 
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Shear stress is also well predicted (see Figure 3) where the standard fc-e model 
over-predicts the shear stress at the edge of the jet in comparison to the data of 
Gutman and Wynanski* and Bradbury’. The k-e-S model presets a narrower 
profile, also indicated by the smaller spreading rate, which agrees with Bradbury s 
measurements. Towards the centerline of the jet, both models accurately predic 

the increase in shear stress and the extrema near \ = 0.08. 

For the next two figures, Figures 4 and 5, there is no experimental data available 
for comparison. The dissipation reaches a maximum at the location of the peak 
shear stress, Figure 4. From this extrema, the dissipation decreases toward the 
centerline and also eventually trails off to zero at the edge of the jet. igure 
shows the behavior of the time scale for the planar jet. Notice that the minimum 
energy transfer time between the large scale structure to the dissipative microscale 
is in the middle of the jet (where S is a maximum since the dimension of S is t ). 

2.4.2 Axisymmetric Jet 

A primary problem with the standard k-e model is its predictive capability for 
a 2-D versus a 3-D jet, the free jet anomaly. In Figures 6 - 10, the standard k e 
model is unable to predict any of the turbulence quantities correctly, overpredicting 
not only the turbulent kinetic energy and the shear stress but also the mean velocity. 
The standard k-e can be “fixed” using a suggestion by Pope 9 . Writing the standard 
e equation for reference, 
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the correct spreading rate can be obtained by changing Ci from its standard value of 
1.45 to 1.6 (see Table 2). Clearly, by modifying a coefficient, the standard k e model 
is capable of predicting the correct spreading rate; however, this model cannot be 
considered general. The spreading rate predicted by the k-e - S model is very close 
to the experimentally measured rate and has been obtained with no modifications 

to the model. t in n 12 17 - 

Comparisons between the k-e - S model and the experimental data , igur 

6 - 10, indicates that the new model accurately predicts the mean velocity profiles 
and the shearing stress across the jet. Less accurately predicted by the fc-e-S mo e 
is the centerline turbulent kinetic energy, Figure 7, where there is considerable 
scatter in the experimental data. 


2.4.3 Planar Mixing Layer 

For the case of a planar mixing layer, neither model accurately predicts the high 
speed side of the layer (see Figures 11 - 15). The experimental data contains 
a much more dispersive edge compared to the computations, Figure 11. One ex- 
planation may be wind tunnel noise. Adding a boundary condition specifying a 
given noise level, 12 percent fluctuating velocity, increases the spreading on the 
high speed side of the computations and improves the comparisons with the data. 
This boundary condition is used with both models. Unfortunately, there is no way 
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to determine what the true noise level was in the experiments or whether its effect 
has been properly accounted for computationally. 

An analytical solution for the mean velocity profile has been described by Schlichting 14 . 
This analytic solution lies between the experimental data and the two computational 
models for the case when a = 9.5 {a is a “tweeking” constant used to tune the ana- 
lytic solution to match the experimental data. Typical values of a are between 9 13 
and ll 15 .) 

Both the k-e-S and the standard k-e model predict the peak value of turbulent 
kinetic energy quite well (Figure 12). However, this peak is shifted toward the low 
speed side of the mixing layer compared to the experimental peak. Patel 16 indi- 
cated that the data of Wygnanski and Fielder 13 had shifted towards the low speed 
side, which is opposite from the trend noticed in this study. Since the discrepancy 
between the data and the computations lies on the high speed side of the mixing 
layer, the problem could be a result of wind tunnel noise as previously mentioned. 

In the the plot of Reynolds stresses, Figure 13, both models over predict the shear 
stress by nearly 30 percent and consequently over-predict the spreading rate given 
in Table 2. 

The fc-e- S model and the standard k-e model predict almost identical levels of 
dissipation (Figure 14). Since a turbulent kinetic energy level is specified at the high 
speed boundary, the dissipation is also specified on this boundary. A constant value 
of = 6.5 x 10 -6 gave the best fit to the data. Figure 15 shows the inverse time 
scale* cfistribution through the mixing layer. The minimum energy transfer time, 
where 5 is a maximum, corresponds to the location of the maximum shear stress. 

At the high speed edge of the mixing layer, the time scale is specified to be S = | . 

2.4.4 Planar Wake 

For the case of the planar wake, the mean velocity profile is correctly predicted 
towards the centerline of the wake, but drops off too quickly at the wake edges (see 
Figure 16). Although the mean velocity shows the correct shape in Figure 16, the 
centerline velocity is high compared to the analytical solution described in Reference 

dy i 

[14] (see Figure 17). As a result, the spreading rates, ( u^-Ul Hz )> predicted by 
the turbulence models are 6 percent too low for the k-e - S model and 3 percent too 
low for the k-e model as compared to the analytical solution. In other words, the 
turbulence models are predicting a much more compact wake with a smaller velocity 
defect than is seen in the experimental data or the analytical fit. This trend is also 
seen in the turbulent kinetic energy profiles in Figure 18 where the k-e - S model can 
predict the correct maximum intensity but drops off too quickly towards the edges. 

Also, although the computed peak value for shear stress is correct the width of the 
curve is too narrow, with the location of the extrema lying too near the centerline 
(see Figure 19). 

The dissipation levels computed by the models apparently reduce the turbulent 
intensity of the flow too much near the edges of the wake, which prevents the wake 
from spreading. There is unfortunately no experimental data to compare with the 
dissipation curves in Figure 20. Interestingly, the location of the minimum energy 
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transfer time in Figure 21 is towards the outer edges of the' wake where the mean 
flow should exert its greatest influence. 



Experiment 

k-e-S Model 

k-e Model 

Planar Jet 

0.11-0.12 

0.114 

0.114 

Round Jet 

0.085-0.095 

0.095 

0.126, ci = 1.45 
0.0928, ci = 1.6 

Planar Mixing Layer 0.16 

0.186 

0.186 

Planar Wake 

.101 14 

0.095 

0.098 


Table 2. Spreading Rate Comparisons for Free Shear Flow 


3. Future Work 

Future plans include: _ „ 

i) Extending this model to wall bounded flow. Specifically, a low Reynolds number 

version of this model will be developed for near wall turbulence. 

ii) Continuing to run this model for more test cases over a variety of flow fields. 


4. 

1 

2 
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Figure 1: Mean velocity profile for a turbulent, planar 
jet. U m az'- centerline velocity. 



Figure 3: Shear stress profile for a planar jet. Umax'- 
centerline velocity. 
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Figure 5: Time scale profile for a planar jet Figure 7: Turbulent kinetic energy profile for an ax- 

isymmetric jet. Umax' centerline velocity. 



Figure 6: Mean velocity profile for a turbulent, axisym- pjgu re g ; Shear stress profile for an axisymmetric jet. 
metric jet. U ma x ■ centerline velocity. U mat : centerline velocity. 





Figure 11: Mean velocity profile for a planar mixing 
Figure 9:. Dissipation rate profile for an axisymmetric iayer u mail centerline velocity. — analytical solution 
jet. Umax- centerline velocity. 



Figure 12: Turbulent kinetic energy profile for a planar 
Figure 10: Time scale profile for an axisymmetric jet u^ng layer. Umax- centerline velocity. 
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Figure 15: Time scale profile for a planar mixing layer 

Figure 13: Shear stress profile for a planar mixing layer. 

U m ax : centerline velocity. 






Figure 17: Axial mean centerline velocity for a planar 
wake behind a flat plate. U m : freestream velocity. 


Figure 19: Shear stress profile for a planar wake be- 
hind a flat plate. Uctuur- centerline velocity, 
freestream velocity. 
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